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This work presents accurate ab initio determination of the hyperfine splitting for the ground 

state and a few low-lying excited states of ^^^Cd+; important candidate for the frequency 

standard in the microwave region, using coupled-cluster theory (CC) in the relativistic frame- 

^ ' work. The hyperfine energy splittings, calculated first time in literature, are in good agree- 

2 ' ment with the recent experimental results. We have also carried out the lifetimes of the 

^ , 5p^Pi/2 and 5p^P^/2 states, which are in good agreement with the available experimcn- 

i tal results. The role of different electron correlation effects in the determination of these 

C/3 ■ 

quantities are discussed and their contributions are presente in the CC terms. 

^ : 

I. INTRODUCTION 

(N ■ 

> : 

I The current frequency standard is based on the ground state hyperfine transition in ^^^Cs 

■ which is in the microwave regime and has an uncertainty of one part in 10^^ Trapped and 
I laser cooled ions are excellent candidates for many high precision measurements Due to 

■ the decoupling of the internal states caused by the perturbations arising from the collisions and 
^ ! Doppler shifts, trapped and laser-cooled ions have been regarded as nearly isolated quantum 

^ ' systems |4|. Therefore, precise measurement of transition frequencies of several species of 

trapped ions have been performed for the purpose of developing the better frequency standard in 
microwave and optical frequency regions J, l5|]. Singly ionized cadmium (Cd"'') has a potential 
for the applications in quantum information processing where the microwave transition between 
the hyperfine states of the ground state is used for both a re-pumping process and manipulation 



of quantum states of trapped ions Recently, ^^"^Cd^ has been proposed for the design of 
a space , a series of measurements for the ground 



state hyperfine splitting of ^^^Cd"*" have been performed 



In recent experiments, the 



control of quantum states of ^^^Cd^ with high degree has become possible [2i] because it has 



2 



a simple energy- level structure and an accessible wavelengths for excitations. This can also 
be regarded as a step towards the development of a new frequency standard in the microwave 
regime Jj. In this work, we have used the coupled-cluster theory with single, double and 
partial triple excitations (CCSD(T)) in the relativistic frame work to calculate the hyperfine 
splitting of ground state and few low-lying excited states of ^'^^Cd'*'. This is the first ab-initio rela- 
tivistic many-body study of the hyperfine splitting of ^^^Cd"*" in literature, known to our knowledge. 



The hyperfine interaction in an atom is generated due to the interaction of different electromag- 
tic 1 
by B 



netic multipole moments of the nucleus and of an atom. The corresponding Hamiltonian is given 



i//,/, =^M«.tW, (1.1) 

k 

where M^") and T^^) are the spherical tensor operators of rank k of the nucleus and the elec- 
tronic system, respectively. Since, nuclear spin of ^^^Cd"*" is ^, hyperfine splitting due to electric 
quadrupole moment will be zero. Magnetic hyerfine interaction, which is the only important hy- 
perfine interaction in this system, for a relativistic electron with nuclear magnetic moment iij is 
given by 



Hhfs = 2^ eca; ■ ^3 (1.2) 



a\ are Dirac matrix for the ith electron. 

In the first-order perturbation theory, the hyperfine energies E^fsiJ) of the fine-structure state 
\JMj) are the expectation values of the hyperfine interaction Hamiltonian. Details of the expres- 



sion are given in Cheng and Childs [s]. The energies corresponding to magnetic dipole hyperfine 
transitions are defined as 

Emi = AK/2 (1.3) 

Here K = 2(7 • J) = F{F -|- 1) — /(/ + 1) — J( J -|- 1) with / and J are the total angular momentum 
of the nucleus and the electronic state, respectively. F is the total angular momentum of an atom 
(nucleus -|- electron). The hyperfine constants A corresponding to the magnetic dipole hyperfine 
interactions are given by 

1 (JWT^'^MlJ) 

A=-iJiH,MJ)=,..,^jyy^d=, (1.4) 
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where /Utv is the Bohr magneton, gj = fij/I where fij is the nuclear dipole moment. The T^^^ 
operator is defined as 

T(i) = J2 -le^^rfa, ■ Yf\n) (1.5) 

i 

with y['^'* as vector spherical harmonic. 

The Fock-space Multi-reference Coupled Cluster (FSMRCC) theory for one electron attachment 
process used here has been described elsewhere [ll|, 12, 13|. We provide a brief review of this 
method. The theory for a single valence system is based on the concept of common vacuum for 
both the closed shell N and open shell Nil electron systems, which allows us to formulate a direct 
method for energy differences. Also, the holes and particles are defined with respect to the common 
vacuum for both the electron systems. Model space of a (n,m) Fock-space contains determinants 
with n holes and m particles distributed within a set of what are termed as active orbitals. For 
example, in this present article, we are dealing with (0,1) Fock-space which is a complete model 
space (CMS) by construction and is given by 



K'^)) = EC..I#^^) (1.6) 

i 

where Cj^'s are the coefficients of and (^^'^'^^'s are the model space configurations. The 

dynamical electron correlation effects are introduced through the valence-universal wave-operator 



Id, 



m 



n = {exp(5)} (1.7) 

where 

m n 

s = j2T. s^'''^ = + + s^'''^ + ■■■ (1-8) 

k=0 1=0 

At this juncture, it is convenient to single out the core-cluster amplitudes S^^'^^ and call them T. 
The rest of the cluster amplitudes will henceforth be called S. Since is in normal order, we can 
rewrite Eq. (|1.7p as 

n = exp(r){exp(5)} (1.9) 

In this work, single (Ti, 5i) and double excitations {T2, S2) are considered for T and S clusters 
operator. Wavefunction of the system with single valence orbital v 

= n^^DP) = e^i+^Hl + + S2v}\^df). (1.10) 
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Triple excitations are included in the open shell CC amplitude which correspond to the correla- 
tion to the valence orbitals, by an approximation that is similar in spirit to CCSD(T) [l^. The 
approximate valence triple excitation amplitude is given by 

"■^'^ ~ r -Lc -Lc _c -c- -r ' 
ta \ \ 

where S'^^'^^^H. are the amplitudes corresponding to the simultaneous excitation of orbitals a, b, k to 
p,q,r, respectively; VT2 and are the connected composites involving V and T, and V and 

^(0,1)^ respectively, where V is the two electron Coulomb integral and e's are the orbital energies. 
The expectation value of any operator O can be expressed, in the CC method, as 

{^v\0\^v) 



^ {<^,\{1 + Sj}e^^0e'^{l + S,}\^,) 
{^,\{l + Sj}eT^eT{l + S,}\<^,) 

The contribution from the normalization factor is given by. 



(1.12) 



Norm = {^^\0\^^){^ -1} (1.13) 

with = (<I?„|e'^^e'^ + {Sle^'' Sv}\^v) for the valence electron v. 

We calculate the DF wavefunctions \^df) using the Gaussian- type orbitals (GTO) as given in 
[isl using basis function of the form 



= Ci^'r'e-'^'''' (1.14) 

with /c = 0, 1, 2, 3, .... for s, p, d, f type orbital symmetry respectively [3]. The radial functions 

'F^' and 'F-^' represent the basis functions correspond to large and small components of the Dirac 

L IS 

orbitals. are the normalization constant which depend on the exponents. The universal even 

tempering condition has been applied to the exponents ; i.e., for each symmetry exponents are 
assigned as 

ai = aQ(3'-^ i = l,2, (1.15) 

where N is the number of basis functions for the specific symmetry. In this calculation, we have 
used ao = 0.00525 and /3 = 2.73. The number of basis functions used in the present calculation is 
32, 32, 30, 25, 20, 20 for / = 0, 1, 2, 3, 4 symmetries, respectively. 
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Number of DF orbitals for different symmetries used in the CC calculations are based on con- 
vergent criteria of core correlation energy for which it satisfies numerical completeness. There are 
only 10, 9, 8, 7 and 5 active orbitals including all core electrons are considered in the CCSD(T) 
calculations for / = 0, 1, 2, 3, 4 symmetries, respectively. We first calculate T amplitudes using the 
CC equations of closed shell systems and then solve the S amplitudes from the open shell equation 
for this single valence states of Cd^. 

TABLE I: Radiative lifetime(ns) for different low-lying states of ^-'^•^Cd+. 



State 



Experiment 



Other theories 



This work 



1/2 



3/2 



3.2(2),'^ 3.05(13),'' 3.11(4),'= 
3.5(2),^= 4.8(1.0),'' 3.14(0.011)^= 
2.5(3),^^ 2.70(25),'' 2.77(7),^= 
3.5(2),^= 3.4(7),'' 2.647(0.01)'= 3.0(2),^ 



2.92,'^ 2.99,'^ 2.74,'' 
3.11^ 
2.50," 2.3,'' 
2.77^ 



3.093 



2.602 



"^Time-resolved laser-induced fluorescence 
''Hanle-theory 



1^ 



"^Beam-laser, beam foil (ANDC) 
"^Phase-shift method 
*=Ultrafast laser pulses 



^Many- 
9Hanle 



)ody third order perturbation theory 



23| 



22 



We report our lifetime results along with the other calculated and measured results in table I. 
As seen, there are large disagreements among earlier estimations. The most reliable experimental 



results are presented by Moehring et al. 2l| to date with a total uncertainty of 0.4%. We can 
see from table I, our calculated lifetime for 5P fine structure states of ^^■^Cd^ are in excellent 
agreement with the recent measured result 21[ and highly correlated MBPT calculations 22l |. 



Table II presents the computed values of the hyperfine energy splittings for ground state and 
few low-lying excited states of ^^^Cd^. We have used the expression (1.3) to compute the highly 
sensitive property to the electronic wavefunction near to the nuclear region. To calculate the 
hyperfine splitting constants corresponding to the magnetic dipole A with /x/ = -0.8278 2j]. It is 



clear from the table II that our calculated hyperfine energy splittings are in excellent agreement 
with the measured hyperfine energy splittings wherever available. We have also estimated the 
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hyperfine splitting of few other excited states for which there are no experimental results available 
to the best of our knowledge. 

All the core orbitals are considered as active in our calculations. In Table III, the individual 
contribution from the one body and effective two body terms of the magnetic dipole hyperfine 
structure constant A for ^^^Cd^ are listed. The first term (O) is the Dirac-Fock (DF) contribution. 
From the differences of DF and total CC results, it is evident that, the electron correlation effects 
to the calculated A value vary from (15-45)% among different low-lying states. 

We know, the Brueckner pair-correlation effects are in the form of O Siv and its conjugate terms 
whereas core-polarization effects are in 52,; and its conjugate terms, in its lowest order. Both 
the correlations are important in the precise determination of the final results. Prom table III, it is 
clear that the largest contribution of electron correlation to the hyperfine splitting constants A for 
different low-lying states comes from the pair-correlation effects. However, the corc-polarization 
contributions are not that least significant. In figure 1, we have plotted the important correlation 
effects to the hyperfine constants of different states with respect to their DF contributions. The 
core-correlation effects seem to be relatively small compared to the other two effects presented in 
the figure. The percentage contribution of the former one is almost the same for all the considered 
states. We see the ratio of pair-correlation and core-polarisation effects is almost one for 651/2 
state, which is different for other states. The pair-correlation effects for the fine structure state 
of 5P arc extremely strong (almost 25%) with respect to DF value, whereas for other states it is 
around 10%. 

Among the other correlation effects, the prominent contributions are observed from S^^ O S2v + cc, 
which is almost 2% for the S and P1/2 states, and more than 5% for P3/2 states. Contributions 
from the effective two-body terms are also significant and comparable with some of the one body 
effects like O S2v + cc for all the states. 

II. CONCLUSION 

In this work, we have determined the hyperfine energy splittings of the ground and a few low- 
lying excited states of ^'^'^Cd'^ using the coupled-cluster theory in relativistic framework. Lifetimes 
of 5P fine structure states of system are estimated also. Our results are in excellent agreement with 
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TABLE II: Hypcrfine energy splitting of different low-lying states of ^^'^Cd 



State Experimental This work 

551/2 15.2(2 Hz) GHz [4, 6] 15.28 GHz 

5Pi/2 2.45 GHz [4] 2.43 GHz 

5P3/2 800 MHz [4] 812.04 MHz 

6S'i/2 3.23 GHz 

6P1/2 667.81 MHz 

6P3/2 236.24 MHz 



TABLE HI: Contributions of different coupled-cluster terms to the ^^■^Cd+ magnetic dipole {A) hypcrfine 
constant. CC stands for the complex conjugate part of the corresponding terms. 

Terms 5si/2 5pi/2 5p3/2 6si/2 6pi/2 6P3/2 

state state state state state state 






-11986.36 


-1837.36 


-284.56 


-2753.76 


-553.68 


-87.78 





-11896.12 


-1823.30 


-284.76 


-2734.74 


-550.04 


-87.82 


Slv + CC 


-2266.42 


-476.44 


-74.19 


-248.89 


-75.64 


-12.52 


S2v + CC 


-1067.66 


-107.83 


-27.36 


-223.06 


-38.22 


-9.47 




-107.95 


-31.23 


-4.85 


-5.66 


-2.65 


-0.45 


S\y S2V 4 


' CC -83.28 


-16.14 


-3.95 


-3.99 


-3.16 


-0.81 


^2v S2V + 


CC -338.50 


-45.99 


-21.23 


-81.69 


-11.16 


-9.13 



Important effective two-body terms of O 





h CC -65.28 


-9.57 


-1.43 


-15.64 


-2.82 


-0.43 


SLOT2' 


h CC 156.55 


23.20 


3.25 


33.77 


5.79 


0.80 


Norm. 


374.86 


51.78 


8.36 


42.14 


10.00 


1.73 


Total 


-15285.99 


-2434.01 


-406.02 


-3237.10 


-667.81 


-118.12 



the available measurements. This suggests the robustness of the CC method and our numerical 
approach in obtaining the accurate wavefunctions of the system considered. 
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5s'l/2 5p'l/3 5p'3/2 ' 6s'l/Z ' 6p'l/Z ' 6P '.3/2 

Atomic States 



FIG. 1: The ratio of core-correlation, pair-correlation and core-polarization effects w.r.t. the DF values 
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